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PROJECT  DESCRIPTION 


I 

I  This  project  concerns  the  demonstration  of  the  leasibility 

I  of  the  tuning  and  stimulation  of  nuclear  radiation.  It 

represents  a  critical  line  of  investigation  in  our  overall 
program  concerned  with  the  feasibility  of  a  gamma-ray  laser. 

Theory,  supported  by  our  experiments  conducted  uncer  this 
contract,  has  indicated  that  anti-Stokes  Raman  upconversion  of 
intense  but  conventional  laser  radiation  proauced  by  scattering 
from  isomeric  states  of  nuclear  excitation  could  lead  to 
significant  sources  of  tunable  gamma  radiation  characterized  by 
the  natural  Mossbauer  widths  of  the  lines.  .^urtner  computations 
have  suggested  that  this  type  of  coherent,  as  well  as  a  type  of 
incoherent,  optical  pumping  could  even  lead  to  appreciable  levels 
of  inversion  of  the  populations  of  nuclear  levels.  thus 
supporting  the  growth  of  stimulated  gamma-ray  intensities. 
Whether  or  not  these  processes  can  reach  thresnold  depends  upon 
tne  resolution  of  basic  issues  that  nave  not  oeen  previously 
addressed  in  an  interai scipi inary  region  oetween  guantum 
electronics  and  nuclear  pnysics.  It  is  the  purpose  of  this 
contract  work  to  study  these  issues  experimentally  in  order  to 
guide  the  development  of  the  technology  and  methods  neeaea  to 
exploit  the  enormous  potential  of  t.his  effect. 

SCIENTIFIC  PROBLEM 

The  viability  of  the  concept  for  tne  tuning  of  gamma 
radiation  by  adding  the  variaoie  energy  or  an  optical  photon 
produced  by  a  tunable  laser  aepencs  upon  the  existence  in  the 
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nucleus  of  a  particular  arrangement  oi  excited  states.  A 
suitable  energy  difference  would  would  make  it  possible  to  dress 
the  nuclear  states  with  the  laser  photons.  Transitions  between 
the  dressed  states  would  then  occur  at  tne  sum  and  aiiference 
frequencies  characteristic  of  the  nuclear  transition.  plus  or 
minus  the  energies  of  integral  numbers  of  laser  pnotons. 

Whether  the  necessary  arrangements  ot  nuclear  states  do 
exist  is  the  central  issue  being  adaressec  in  this  contractea 
work.  Surprisingly.  such  inlormation  is  currently  unknown 
because  such  potentially  useful  states  would  lie  in  tne  blinc 
spots"  ot  the  conventional  techniques  ot  nuclear  spectroscopy. 
Normal  Mossbauer  spectroscopy  provides  enormous  resolution,  out  a 
tuning  range  that  is  inadequate  by  orders  of  magnitude  to  support 
any  possible  study  of  transitions  to  the  intermediate  states  ot  a 
multiphoton  process.  Conversely,  crystal  spectrometers  provide 
broad  tuning  ranges.  but  levels  of  resolution  that  miss  py  two 
orders  of  magnitude  the  threshold  tnat  would  be  necessary  to 
separate  the  transitions  to  the  initial  and  intermediate  states. 
As  a  consequence.  the  ideal  arrangement  ot  nuclear  energy  levels 
needed  for  the  Raman  upconversion  process  couia  be  a  common 
occurrence  that  has  gone  unnoticed  because  ot  tne  inadequacies  ot 
conventional  nuclear  spectroscopy. 

The  critical  problem  in  this  research  has  two  tacets.-  l; 
the  development  of  an  appropriate  spectroscopic  technique,  and  2) 
the  search  for  a  suitable  medium  for  a  large-scale  eltect.  The 
dressing  of  the  nuclear  states  not  only  affects  their  energies, 
but  also  changes  their  transition  properties.  .'^orbidden  nuclear 
transitions  should  become  allowed  so  that  the  me t as t ao i 1 i t v  oi 
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isomeric  states  would  be  "switched  oft"  as  the  states  were 
dressed.  This  would  greatly  enhance  the  prospects  :or 
stimulating  the  gamma-ray  transition,  in  addition  to  rendering  it 
tunable.  It  is  the  development  of  the  investigative 
instrumentation  and  the  verification  ot  these  predicted  eftects 
that  comprise  the  scientific  problem  addressed  by  this  contract 
research . 


TECHNICAL  APPROACH 

For  the  resolution  of  the  cent.^al  issue  of  the  existence  of 


potentially  useful  intermediate  states  in  a  multiphoton  upconver- 
sion  of  optical  photons  to  gamma-ray  e.nergies,  it;  was  firs: 
intended  to  demonstrate  sum  frequency  generation  in  one  case  in 
which  nonresonant  intermediate  states  were  known  to  exist.  This 
was  the  case  in  which  both  initial  ana  intermediate  states  were 
magnetic  sublevels  of  the  same  nucleonic  state  and  in  which  tne 
transitions  were  mediated  by  the  Ml,  magnetic  dipole  operator. 
Experimental  data  reproduced  in  the  literature  suggested  that 
such  a  process  had  already  been  unknowingly  demonstrated  for  the 
generation  of  r ad io f r equency  sidebands  to  Mossbauer  transitions 
at  the  sum  and  difference  frequencies.  This  suggested  the  cevei 
opment  of  a  new  instrument,  a  ."requency  .'Modulation  Spectrometer 
for  gamma-ray  energies,  designed  to  support  the  needed  studies  of 
nuclear  structure  with  the  precision  of  .Mossbauer  spectroscopy 
applied  over  a  tuning  range  of  energies  lying  considerably  bevonc 
the  state-of-the-art  at  the  time  our  work  began.  Pursuant  to 
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field  with  measurement  and  parameterization  of  the  conversion 
efficiency  into  the  sum  frequency  line  to  determine  the  practical 
limits  on  the  ultimate  linewidths  and  tuning  ranges  that  can  be 
achieved.  This  technique  will  then  be  used  in  a  "bootstrap" 
approach  to  support  a  search  for  accidentally  resonant  inter.med- 
iate  states.  By  replacing  the  radiofrequency  excitation  with 
tunable  higher  frequencies,  it  is  expected  that  the  tuning  range 
of  Mossbauer  spectroscopy  can  be  extended  by  further  orders-of- 
magni tude . 


PROGRESS  DURING  THIS  REPORTING  PERIOD 
Instrumentation  of  FMS 

The  sum  and  difference  frequency  sidebands  proaucec  on 

intrinsic  Mossbauer  transitions  has  made  possible  very  eifective 

new  instrumentation  for  high  resolution  spectroscopy  at  gamma-ray 

energies.  A  prototype  version  of  this  .'frequency  Modulation 

1 

Spectrometer  (FMSj  was  first  described*  by  our  iaooratory  in 
1985,  and  subsequent  refinements  were  made  during  the  current 
reporting  period  under  this  contract  and  under  a  related  SDI 
contract  directed  by  ONR.  This  device  monitors  changes  in  the 
intensity  of  transmitted  single- f requency  gamma  photons  as  a 
function  of  frequency  of  the  long  wavelength  photons  o:  tne 
alternating  magnetic  field  in  which  the  absorbing  nuclei  are 
immersed . 

During  the  past  year  our  prototype  '.Nuclear  Raman 
Spectrometer"  was  refined  into  more  mature  technology  resulting 
in  the  Frequency  .Modulation  Spectrometer  (FMS)  for  gamma-ray 
energies  (Figure  1).  The  original  prototype  device  nag  required 
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a  tedious  level  of  manual  interaction,  and  this  was  replaced  with 
a  fully-automated  and  computerized  control  system  during  the 
current  reporting  period.  At  its  heart  is  a  multi-channel  scalar 
(MCSl  and  IEEE-488  GPIB  interface  with  an  Apple  11+  computer. 
The  MCS  was  designed  to  have  a  lOOl  duty  cycle.  The  GPIB  enables 
the  spectrometer  to  sweep  continuously  through  the  frequencies  of 
an  rf  magnetic  field  produced  with  a  Wavetek  frequer.ry 
synthesizer.  The  Mossbauer  drive  allows  the  frequency  of  the 
gamma  photon  to  be  biased  by  a  constant  Doppler  shift,  if 
desired.  In  its  present  form,  the  FMS  device  has  an  instrumental 

9 

resolution  of  lOOHz  and  a  continuous  tuning  range  of  10  Hz  with  a 
stability  of  O.lHz/sec  with  no  mechanical  movements  required 
anywhere.  These  characteristics  are  comparable  to  a  Mossbauer 
spectrometer  with  a  means  of  shifting  the  gamma-ray  source, 
having  a  resolution  of  lOnm/sec  and  a  range  of  lOOmm/sec  with  a 
stability  of  0 . 0 Inm/ sec / sec .  Demonstration  spectra  were  acquired 

57 

with  Fe  showing  isomer  shifts  and  t.hermal  shifts.  Because  we 
were  using  a  modulation  type  of  spectroscopy,  the  static  features 
could  be  suppressed,  and  these  different  effects  were  obtained 
with  unprecedented  clarity. 

57 

In  operation,  .'^MS  of  Fe  provides  a  direct  measurement  o: 
rf  sideband  positions  and  intensities  from  which  one  can 
extrapolate  information  about  the  transitions  between  Zeeman 
split  energy  levels  (parent  transitions],  labeled  1  througn  b  in 
Figure  2.  Radiofrequency  sidebands  have  been  labeled  as  a 
parent  transition  preceded  by  a  number  of  +'s  or  -’s,  the  numoer 
of  which  corresponds  to  the  number  of  rf  field  energy  quanta 


(order  of  the  sideband)  added  to  or  subtracted  from  the  parent 

transition.  The  symmetrically  opposed  parent  transitions  1  and  6 

are  separated  by  123.7  MHz  in  pure  iron.  .Applying  a  61.85  MHz 

alternating  magnetic  field  to  the  Fe  foil  produces  L+l}  and  (-6) 

sidebands  which  overlap  in  the  symmetric  center,  or  transition 
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center,  of  the  hyperfine  structure  of  the  Fe.  The  energies  of 
the  gamma  ray  emitted  by  the  source  and  the  transition  center  of 
the  absorber  differ  by  the  isomer  shift.  A  (Figure  3a).  In  !^MS 
the  Stokes  sideband  from  parent  transition  6,  (“6)  would  be 
detected  at  a  frequency  of  (61.85-a)  MHz  while  the  anti-Stokes 
sideband  from  parent  transition  1.  (+1)  would  be  detected  at  a 

frequency  of  (61.85+a)  MHz  (Figure  3b  and  c).  Therefore,  .^M3 
should  produce  a  spectrum  with  two  peaks  around  60  .MHz,  separated 
by  2a  . 

If  we  apply  a  small  Doppler  shift.  3.  to  the  source,  we 
should  obtain  an  FMS  spectrum  with  two  peaks  around  60  MHz. 
separated  by  2(a  d)  (Figure  4).  Classically,  tne  irequencies 
at  which  sidebands  appear,  fs,  is  simply 

fs(MHz)  =  (  V  -  ( P  .  isojJ  ♦  ill.b/ord),  (,1; 

where  P.  is  the  position  of  the  i ’ th  parent  transition  imm/secj, 
V  is  the  velocity  of  tne  source  imm/sec),  iso  is  the  isomer 
shift  (mm/sec),  and  ord  is  the  order  of  the  sideoand  of  interest. 
The  source  used  was  in  a  Pd  lattice  (iso  =  -0.135). 

Rad io f r equency  sideband  positions  are  also  apparentiv 
affected  by  the  intensity  of  the  rf  magnetic  field  (Figure  5). 
It  is  yet  to  be  determined  whether  sideband  position  is  a 
function  of  intensity  as  well  as  frequency  of  the  rf  magnetic 


field,  or  whether  the  temperature  snift  of  the  parent  transitions 


is  being  detected,  or  both.  Since  we  start  with  a  negative 
isomer  shift.  raising  the  temperature  of  the  absorber  snouid 
reduce  the  energy  difference  between  t.de  source  transition  ana 
the  transition  center  of  the  absorber.  Therefore,  increasing  tne 
rf  field  intensity  should  raise  the  temperature  of  the  absorber 
and  in  turn  decrease  the  isomer  shift. 

It  seems  clear  that,  in  addition  to  providing  information 
about  sideband  intensity  and  position,  FMS  could  also  prove  to  be 
a  means  for  direct  and  accurate  measurement  of  isomer  ana 
temperature  shifts,  spectroscopic  quantities  that  are  difficult 
to  measure  with  .Mossbauer  spectroscopy  as  usually  practiced 
because  of  the  difficulty  in  obtaining  such  sma 11  velocities  with 
such  precise  control.  .'^or  the  purposes  of  the  gamma-ray  laser 
program,  it  is  the  combination  of  narrow  instrumental  width  and 
large  tuning  range  that  offers  the  greatest  attractions. 

To  dress  an  isomeric  state  requires  a  certain  arrangement 
of  nuclear  levels  that  would  make  them  undetectable  to 
conventional  techniques  of  nuclear  spectroscopy.  Our  method  ot 

is  the  only  means  found  to  dace  that  can  oe  used  to  sea.-cn 
for  this  comoination  among  the  29  best  candidates  for  a  gam.ma-riv 
laser.  The  successes  ot  t.he  new  FMS  apparatus  for  nuciea.' 
spectroscopy  indicate  that  a  much  nigher  resolution,  ov  per  naps 
several  more  orders  of  magnitude,  can  be  achieved  tnrougn  a 
reasonable  upgrade  of  the  apparatus.  Ii  the  range  of  tunaoilit'.- 
does  extend  to  the  ferromagnetic  spin  resonance  i  FSP. )  frequency, 
will  pe  possible  to  construct  a  swept  frequency  oevic,- 


then  it 
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capable  c:  continuously  tuning  over  a  range  of  10  linewidths, 
an  enormous  improvement  in  the  s t a t e'o f ■ t he • ar t  of  nuclear 
spectroscopy 
Theory  oi  FMS 

Most  .aboratory  sources  ci  gamma  raoiation  emit  at  levels  of 
intensity  corresponding  to  single  photon  conditions.  Mossbauer 
experiments  are  rarely  conducted  at  such  great  intensities  that 
the  detection  of  two  pnotons  would  be  probable  in  the  transit 
time  spent  between  source  and  absorber.  Under  those  conditions, 
the  perception  of  gamma  rays  as  streams  of  particles  is 
instinctive.  Nevertheless.  as  elements  of  electromagnetic 
radiation  they  must  also  be  considered  as  carrier  waves  of  high 
frequency.  As  the  name  implies.  the  Frequency  Modulation 
Spectrometer  derives  its  operation  from  the  modulation  of  the 
gamma-ray  carrier  and  represents  a  translation  into  the  nuclear 
domain  of  one  of  the  powerful  techniques  of  laser  spectroscopy  at 
the  atomic  level.  As  described  above,  the  device  itself  works 
better  than  the  theory  describing  it.  and  during  the  current 
reporting  period  emphasis  developed  upon  refining  and  validating 
the  dressed  state  theory  describing  the  origin  of  the  carrier 
modulation.  Reviewed  in  the  following  material  will  be  five 
different  approaches  to  a  comprehensive  theoretical  description. 
Not  published  elsewhere,  each  is  being  reported  here  as  each 
approach  is  useful  in  defining  acceptable  levels  of 


approxima  t ion . 


CONTEXT  OF  THEORETICAL  EFFORTS 


In  our  version  of  FMS,  the  spectrum  of  the  gamma  photons 

emitted  by  the  source  is  modulated  by  its  passage  through  a  thin, 

2 

ferromagnetic  foil  immersed  in  a  radiofrequency  magnetic  field  . 

57 

Most  work  has  been  done  upon  Fe  with  the  radiofrequency  mag¬ 
netic  field  being  applied  in  the  plane  of  the  foil.  Although  the 
interaction  energy  of  the  nuclear  magnetic  moment  directly  with 
the  applied  magnetic  field  is  insufficient  to  cause  a  significant 
perturbation  of  the  system,  the  applied  magnetic  field  is  suffi¬ 
cient  to  cause  the  direction  of  the  magnetization  within  the 
material,  M(t),  to  process  in  a  complicated  manner  and  this 
interaction  of  the  nuclear  moment  with  M  is  of  significant  magni¬ 
tude  . 

Ferromagnetism  is  a  cooperative  magnetic  behavior  resulting 
from  the  exchange  interaction  between  the  electrons.  In  ferro¬ 
magnetic  materials,  the  principal  effect  of  the  exchange  inter¬ 
action  is  to  cause  parallel  ordering  or  alignment  of  the  atomic 
spins  in  a  common  direction,  so  that  there  is  a  large  spontaneous 
magnetization  even  in  the  absence  of  an  applied  field.  Maximum 
ordering  is  obtained  at  0°K  where  the  randomizing  effect  of 
thermal  agitation  disappears.  At  the  Curie  temperature  the 

magnetic  ordering  is  destroyed  by  thermal  agitation,  and  the 
spontaneous  magnetization  disappears.  Between  these  two  extremes 
a  sample  of  ferromagnetic  material  is  usually  divided  into  small 
volumes  called  domains,  which  may  vary  in  size,  shape  and 
direction  of  magnetization.  Within  each  domain,  the  magnetiza¬ 
tion  is  uniform  and  has  the  maximum  or  saturation  value,  M  , 

s 
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characteristic  of  the  material  and  the  temperature.  Within  a 

ferromagnetic  material ,  the  magnetization  vector  Ml  t )  may  vary  in 

direction  but  not  in  magnitude.  even  when  traversing  a  domain 

wall.  Therefore,  changes  of  the  magnetization  are  precessions 

and  must  be  described  by  equations  which  reflect  this. 

57 

The  magnetic  field  at  the  site  of  the  Fe  nucleus  is  the 

sum  of  several  large  terms  which  may  be  positive  or  negative.  In 

57 

an  Fe  foil  the  magnetic  field  at  the  nucleus  has  been  deter- 

3 

mined  to  be  about  330  x  10  Oersteds  in  the  direction  opposite  to 
the  local  direction  of  the  bulk  magnetization,  M(t).  The  elec* 
trons  respond  to  the  changes  of  M(tj  quite  rapidly  since  the 

magnetic  moments  of  the  electrons  are  strongly  coupled  to  MCt) 
C the  cooperative  effect).  Thus,  in  response  to  any  applied 
fields,  the  magnetic  field  at  the  site  of  the  nucleus  also  pre- 
cesses  in  a  manner  which  reflects  the  precession  of  the  bulk 
magnetization.  For  this  reason,  the  motion  of  the  magnetization 
M(t)  must  be  determined  as  a  function  of  time  as  a  part  of  any 
comprehensive  description  of  the  ultimate  effects  of  applied  .H 
fields  upon  nuclear  states. 

The  equation  of  motion  of  the  magnetization  vector  .*1lt,'  is 
determined  by  the  Landau-Li f shi tz  equation: 

-Cl/|-y()  dMtt)/dt  =  i.M  X  F)  -  A  [F  -  'M-Fl  :  .1  : 

The  gyromagnetic  ratio  of  the  electron  is  given  bv  and  a  is 

an  experimental  parameter  obtained  from  ferromagnetic  resonance 

-  3  - 

studies,  magnetic  field  F  in  equation  3)  is 

the  vector  sum  of  four  contributions:  1  tne  magnetic  fie.o  re- 

t  0 


quired  to  simulate  the  magnetic  anisotropv  effects  inherent 
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the  crystal  structure  of  the  material,  2j  the  "magnetostatic" 

fields  due  to  the  finite  size  of  the  sample,  3)  the  exchange 

field,  which  represents  a  potential  energy  term  leading  to  the 

ordering  of  neighboring  magnetic  dipole  moments  on  a  microscopic 

3  4 

scale,  and  4]  the  applied  magnetic  field.  The  solutions  ’  to 
equation  C2]  depend  strongly  on  the  anisotropy  of  the  material 
and  on  the  geometry  of  the  sample  and  are  usually  of  an  exceed¬ 
ingly  complicated  nature.  Approximations  were  developed  to  vary¬ 
ing  degrees  as  part  of  each  of  the  different  approaches  to  the 
nuclear  sideband  problem. 

During  the  past  year  our  efforts  to  develop  a  model  of  the 
modulation  of  the  gamma  spectrum  have  focused  upon  the  following 
five  approaches,  to  be  summarized  in  this  report  below: 

1.  The  refinement  of  the  theoretical  approach  initiated  by 
C.  B.  Collins  and  B.  D.  DePaola“  in  Optics  Letters  10,  25  C19B5) 
culminating  in  the  results  presented  at  the  1986  IQEC  meeting. 
We  have  referred  to  this  as  the  multiphoton  processes  in  nuclear 
states. 

2.  A  literature  seared  and  study  of  the  applicability  of 
the  'dressed  atom"  approach  of  C.  Cohen -Tannoud ; i . 

3.  The  development  of  a  coherent  switch  model  of  the 
magnetization  M( t )  by  integrating  the  Landau-Li f shi tz  equation. 
This  model  was  specific  for  a  thin  ferromagnetic  foil.  The 
solution  for  .Ml' t  )  is  in  agreement  with  results  presented 
e 1 sewher e . ^ 

4.  The  development  of  a  domain  wall  passage  model.  This 
research  does  not  appear  to  have  been  fruitful,  although  it  may 


still  be  of  use. 


5.  In  collaboration  with  Professor  P.  Berman  of  New  York 
University,  the  development  of  solutions  to  the  time  dependent 
Schroedinger  equation  utilizing  information  about  the 
magnetization  derived  from  approach  3.  The  resulting  method 
appears  to  be  the  best  approach.  The  computer  coae  is  in  the 
process  of  development  and  is  being  verified  first  for  the  static 
case . 

57 

1 )  Multiphoton  Processes  in  Nuclear  States  of  Fe  in  a 
Ferromagnetic  Foil  Immersed  in  a  Radiof reauencv  Magnetic 
ELfl-d 

57 

The  time-dependent  Schroedinger  equation  of  the  Fe  nucleus 
embedded  in  a  ferromagnetic  foil  immersed  in  an  rf  magnetic  field 
is 

=  -  {H^ppCt)-u  M(t)-u  +  H^Ctl-M  ^ 

the  eigenstate  of  the  system, 

the  atomic  Hamiltonian  associated  with  quantum 
number  (  s  J  <x, 

the  magnetic  moment  of  the  nucleus, 

the  applied  radiofrequency  magnetic  field  fa  few 
gauss!  with  frequency  , 

the  magnetization  of  the  material  fkilogauss), 
the  demagnetizing  field  due  to  the  surface  poles 
developed  when  M  rotates  out  of  the  plane  of  the  foil 
Cki logauss  J , 

the  demagnetizing  field  due  to  the  poles  developed 
when  the  nuclear  magnetic  moment  u  rotates  out  of  the 
plane  of  the  foil  (very  small). 


ihCd?/dt) 

where 

H.  = 


H 


app 


M 

H 


H 


dn 
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The  expression  in  the  brackets  represents  the  interaction  oi 
the  nuclear  magnetic  moment,  ju.  with  the  vector  sum  of  ail  of  the 
magnetic  fields  to  which  the  nucleus  is  subjected.  The  interac¬ 
tion  of  the  magnetic  field  with  the  electrons  is  included  in  the 
behavior  of  M  and  therefore  does  not  require  another  term.  This 

notation  and  approach  are  essentially  the  same  as  outlined  by 

2 

Collins  and  DePaola  . 

^  ^  ^ 

Since  the  terms  H  ^  (t)*a  and  H,  (t)‘u  are  small  when  com- 

app  dn 

pared  to  the  other  terms  in  the  Hamiltonian,  they  are  neglected. 

The  simplified  time  dependent  Schroedinger  equation  becomes 
ih  Caf/dt]  =  T  -  {M(t)-il  +  H^Ct)*i!}  f  (4) 

in  the  fixed  laboratory  coordinates. 

A  new  coordinate  system  is  defined  in  the  laboratory,  such 
that  the  magnetization  M  in  the  plane  of  the  foil  corresponds  to 
the  z  axis,  that  is,  as  M  responds  to  the  applied  magnetic  field 
(by  rotating  in  the  plane  of  the  foil),  the  foil  is  rotated  about 
an  axis  perpendicular  to  the  foil  in  such  a  way  that  M  always 
lies  along  the  z  axis  of  the  laboratory  coordinates.  Doppler 
shifts  are  introduced  because  the  foil  is  moving,  but  since  the 
absorption  spectrum  is  observed  along  an  axis  perpendicular  to 
the  foil,  only  the  second  order  transverse  Doppler  shift  occurs, 
which  should  be  negligible.  This  coordinate  system  is  denoted  as 
primed  in  the  following  material.  In  it  M  has  constant  magnitude 
and  direction  (z  axis)  as  shown  in  Figure  6.  Equation  (4)  may  be 
rewr i t  ten 

iM  Of’/dt)  =  f’  -  {M-u  *  H^(t)-y}  If’  (5) 

The  magnetic  moment  of  the  nucleus.  u.  is  related  to  its  angular 


momentum  through  the  equality 
u  - 

Where  x  denotes  either  the  ground  (g)  or  excited  Ce)  state  of  the 
nucleus.  The  symbol  stands  tor  the  nuclear  magneton. 

=  Ceii/2M1  =  0.505  x  10'“^  amp-m^  ,  ',7] 

in  the  MKS  system  of  units. 

The  .Hamiltonian  of  a  particularly  useful  basis  set  of  eigen¬ 
states  and  energies  is  obtained  from, 

Hq  =  H.  -  M*w  =  H.  -  Ww  J  .  I  3  ) 

BA  A  X  2  5 

where  corresponds  to  the  total  splitting  or  the  x  level  by 

the  magnetization  M, 

In  the  basis  set.  the  time  cependent  Schroedinger  equation  of  the 
X  state. 

ih  cay^/ati  =  j,  y;  .  cio] 

has  solutions  of  the  form. 

f^fcx  ,J’,m’3  =  !  cc^ .  J  ’  ,m’ >expf.  -  i  l  /h  -  m’.j  J  t  ]  .  ill) 

X 

The  quantum  numberCs]  describe  the  nuclear  state  x  witn  energy 

E  ,  in  the  absence  of  any  magnetization  riela  in  the  materia:  or 
“x 

any  applied  field.  The  angular  momentum  state  of  the  nucleus  in 
the  primed  coordinate  system  is  described  by  J’  and  ra'.  Tnese 
basis  states  of  the  system  might  be  describea  as  the  ncn-inter- 
acting  basis"  states. 

The  interacting  state  f’  can  De  derivea  from  tne  oasis  state 
by  means  of  a  rotation  througn  an  angle  r  about  an  a.xis 
parallel  to  H^Ct).  This  rotation  operator  is  Hermitian  a no 
unitary  and  of  the  form 
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R  =  exp[i«^n*J]  =  exp(.i«^J^3  .  1123 

where  n  is  a  unit  vector  parallel  to  the  axis  oi  rotation.  The 
interacting  state  is  given  by 

f  ^  =  e'-^*x^d^  ^ 

Combining  equations  (5)  and  (133.  the  time  dependence  of  the 
angle  0^  is  given  by  the  equation 
(;d(^^/dt3  =  9x^N  ^d*-*^^^** 

If  Hj[t3  is  known.  the  angle  can  be  computed  numeri¬ 

cally  and  therefore  may  be  considered  a  known  function. 

The  transition  amplitude  is  determined  bv 


iMCda  /dt3  =  iha 

eg  eg 


115  3 


/r(-*-ii^J.3-.,-  f,  .-(l*,/  ■i'*»tr(-^i(>J_3,-.,  ,,  ..,'1 

/ 1  e  e  d  ,m  31  V^e  y  le  g  d  ,m  3J 


dr 


Since  R  is  Hermitian  and  unitary,  this  equation  may  be  rewritten 


ih  a 


eg 


where 


iw  t  ir  r  ^i  C«^-<<»^3  J  . !  .  1*  ,, 

y/e  e  g  d  flat, J,m3|W 

Jl  o  e  J  ■ 


116  3 

.J’,m'3  dr 


W  =  e'-*^'*g'^d^  V  e''^^*g^d^  V  tl?! 

y  y  V 

The  operator  J.  is  not  diagonal  in  the  primed  coordinate  svstera 
d 

A  ^ 

since  n  is  perpendicular  to  f  the  z  axis;.  It  is  necessary  to 
use  the  finite  rotation  matrix^*^  to  transform  into  a  system  in 
which  the  z  axis  is  parallel  to  n.  After  operating  witn  then 

the  finite  rotation  matrix  is  again  used  to  retur.n  to  tne  primed 
laboratory  system. 

The  time  derivative  of  the  total  transition  amplitude  in 

terms  of  <f’lm  3lU  iV’fm  3>  is  given  bv 
e  e  y  g  g 


ih  a  =  (1/2)^  (cq  exp[in,t]  <-3/2|W  |-l/2>  R,(t) 
eg  D  1  7  1 

Cg  expfin2tl  <  - 1 /2  I  W_^  I  -  1 /2> 

+  Cg  expfiOgt]  <-*-l/2  |W_^  I -l/2>  (18^ 

■*■  expfiO^t]  <- 1/2  !U^  I  ■►l/2>  ^2*-^^ 

■*•  Cj^  expfiOgt]  <  +  1/2  !  I  ■*-l/2>  ^2^’'^ 

+  expfiOgt]  <  +  3/2  I  W_^  I +l/2>  R^lt)) 

where  is  the  detuning  of  the  transition  as  shown  in 

Figure  7 .  That  is , 

n.  =  CaE./K)  -  u  .  [19) 

1  1  > 

where  is  the  energy  of  the  transition.  The  concentra¬ 

tions  of  the  ground  state  levels  are  c^  and  Cg  for  m  =  -1/2  and 
m  =  +1/2  respectively.  If  the  angle  ♦  is  defined  by 


^  9  ^  ^ 

e  g 


(20) 


the  factors  R^[t)  and  ^2^*'^  which  appear  in  the  transition  ampli 


tudes  are  defined  as  follows: 

R^Ct)  =  (1  +\/3j  +  (3  ->/3)  e^ 


R.,(t)  =  [3  +V3 


-iC3«/2) 


+  ( 1  -  y/3  )  e 


i  (^  /  2 


The  transitions  probabilities  are  calculated  from 

rC'iW 

t  ®  t  eg  eg 

When  Rjlt)  and  R2^^)  3’’®  expanded  in  Fourier  series 


(  21a  ) 


[21b) 


R.(t)  =  I  [R  ) 


[23  J 


there  are  two  components  to  the  transition  probability  wnich 


survive  the  time  averaging  procedure.  The  first  is  the  "normal 


contribution  given  by 
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^  c^l<  +  3/2!W^!-^l/2>|  ^  tflOg-nw^jJ 
■"  [  CbK-1/2  IW,.  | -l/2>  !  ^  ilO^-nwJ 


>  c^|<  +  l/2|W^|-l/2>!^  <5in3-nu,i) 


c^|<-l/2(W^|-^l/2>|  ^ 


c^|<  +  l/2|W^|+l/2>l^  tfCOg-nw 


^>]| 


124  ) 


After  R^Ct)  and  ^2 ^ ^ ^  have  been  expanded  in  terms  of  their 
Fourier  components,  a  typical  term  in  the  product  of  the  two 
expansions  shown  in  C22]  is  of  the  form 


c.tME),  C'CME)* 

m ,  n 


^i Cn, -nw) t . , 
e  i. 

^ - i ( n . - mw } t . , 

n .  -  nw 

n .  -raw 

■»  *  • 

]  J 

.  I  2  5  } 


where  (ME)^  is  the  matrix  element  of  the  i''*^  transition.  When 
the  time  average  is  taken  as  shown  in  expression  (22),  this  term 
is  zero  for  all  m  /  n.  The  summation  over  m  collapses  to  one 
term  which  does  not  vanish  when  the  time  average  is  taken.  When 
m  =  n,  the  equality^*^  in  the  limit  t"*®. 


1 im  sin  at 

t  2 

TrtcL 


=  <5  (  a  ) 


(26) 


may  be  used,  which  introduces  the  delta  functions  of  frequency. 

In  reality,  the  delta  functions  of  frequency  must  be  re¬ 
placed  by  Lorentzian  distributions  centered  at  the  position  of 
the  delta  function.  The  width  of  the  Lorentzian  corresponds  to 
the  lifetime  of  the  excited  state.  The  parent  lines  of  the 
transition,  that  is,  the  spectrum  in  the  absence  of  the  rf  raagne- 
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tic  field,  correspond  to  the  n  =  0  amplitudes, 

=  ^(2)’ [c||<-3/2IU^|-1/2>|2  mOj) 

+  c^l<*3/2lW^|-*-l/2>l^ 

+  1  (R2)o1^  5C02)  L27) 

Cg!<-^l/2  IW^I -1/2>|  ^ 

+  c?!<-l/2!W  I+1/2X  ^  <5Cn^  ) 

A  V  4 

+  c^|<  +  1/2|W^|-^1/2>i2  <5105) 

The  delta  functions  of  frequency  indicate  that  the  parent  lines 
will  be  observed  at  their  normal  frequencies, 

0.  =  detuning  =0  or  aE. /h  =  w  .  C28) 

1  1  V 

The  intensity  of  the  first  order  sideband  may  be  found  in 
the  same  manner  by  taking  either  n  =  il.  For  n  =  -^1,  the 
frequencies  of  the  first  order  sidebands  in  the  absorption  spec¬ 
trum  is  found  by  taking 

fl^  -  =  (AEj^/h)  -  -  0  ,  (.29) 

or 


r (ae. /h)  -  w, ]  =  w 

1  1  7 


(30  3 


Thus,  n 


■*■1  corresponds  to  the  first  negative  sideband  and 


conversely  n  =  -1  corresponds  to  the  first  positive  sidebana. 

The  absolute  magnitudes  of  (R. 3  govern  the  amplitudes  of  the 

1  n 

sidebands.  For  this  reason.  it  is  vitally  important  to  have  an 
adequate  model  of  the  behavior  the  magnetization  M  as  a  function 
of  time. 

The  second  component  is  a  quantum  enhancement  wnich  acds  on 


to  the  "normal"  curve  whenever  the  condition. 

0.  -  0.  =  pw,  (313 

1  ]  ^  1 

(p  is  any  positive  or  negative  integer;  is  satisfiec,  that  is. 


>.-  >  .“-i  '  V 


I*-  L’. 


fcV  _%  «  ,.*•  _*w  •  _  •  _‘v  _  ^ 


•  4' 


i>K 
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whenever  two  transition  energies  differ  by  an  integral  number  of 
photons  in  the  rf  field.  The  intensity  given  by  the  "normal" 
curve  is  increased  by  the  amount 


^  {if  ^  '32) 


where , 

K  =  c.  c*  CME).  CME)*  .  133) 

1  ]  1  3 

As  can  be  seen  from  expression  (33),  the  enhancements  to  tne 
absorption  spectrum  fall  on  the  same  frequencies  as  the  usual 
sidebands,  including  the  parent,  and  occur  only  when  equation 
(31)  is  satisfied.  If  i  =  1  or  6,  then  For  all  other 

values  of  the  index  i,  then  index  j  follows  the  same 

rules.  For  this  reason,  the  quantum  enhancement  must  be  deter¬ 
mined  on  a  case-by-case  basis. 

At  this  point  the  model  is  a  direct  extension  of  the  one 

2 

originally  proposed  by  DePaola  and  Collins  which  was  limited  in 

validity  to  the  vicinity  of  a  moving  domain  wall  where  the  angle 

of  rotation  of  the  magnetization  M(t)  in  the  plane  of  tne  foil  is 

large.  The  extension  realized  in  this  section  as  characterized 

in  equations  (24)  and  (32)  is  more  generally  vaiic. 

Unfortunately  both  models  share  the  dependence  on  the  rather 

simplified  description  of  the  magnetization  which  strictly 

confines  M(t)  to  a  plane.  Not  surprisingly,  neither  one 

adequately  describes  the  enhancement  of  the  sideband  intensity 
57 

observed  in  Fe  near  the  applied  frequency  of  45  MHz.  This 
enhancement  of  the  sideband  amplitude  was  thought  to  be  due  to 
the  mixing  of  the  Zeeman  levels  of  the  ground  state  by  the 


applied  magnetic  field.  For  a  better  understanding  of  the  mixing 


effects  of  an  intense  field,  the  ’dressec  state  t.neory  of  C. 

5  6  7  3  9 

Cohen-Tannoud j i  ’  ’  and  S.  Haroche  ’  was  stuaiec.  because  it 
had  been  refined  and  experimentally  verified  to  a  very  nigh 
degree  of  accuracy  for  analogous  cases  in  atomic  physics. 

Also,  before  proceeding  further,  it  is  necessary  to  aeveiop 
a  better  model  for  the  behavior  of  the  magnetization  t)  in  tr.e 
presence  of  an  applied  rf  field.  To  accomplish  this.  it  is 
necessary  to  integrate  the  Lanaau-Li f shi t z  equation  [equation  1] 
as  discussed  in  sections  3  and  4. 
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2 )  “Dressed  State"  Theory 

6  7 

At  first,  the  work  of  C.  Cohen-Tannoud  j  i'' '  ’  and  S. 

8  < 
Haroche  ’  on  atoms  immersed  in  a  strong  electromagnetic  fiela  I 

t 

of  a  single  frequency  (optical  or  radio  frequency)  is  very  | 

j 

attractive  conceptually.  It  departs  considerably  from  tne  usual 
approach  in  which  the  interaction  of  the  atom  with  the  applied 
electromagnetic  field  is  considered  as  a  perturbation  to  be 
handled  using  standard  perturbation  theory.  The  need  for  an 
alternative  arises  from  the  fact  that  perturbation  theory  is  not 
appropriate  when  the  strength  of  the  interaction  is  large  due  to 
the  presence  of  very  strong  fields.  In  the  dressed  state  theory. 

I 

the  quantized  system  consists  of  both  the  atom  and  the  electro¬ 
magnetic  field.  The  electromagnetic  field  is  handled  as  an 
integral  part  of  the  system,  not  a  perturbation,  and  it  is  con¬ 
sidered  quantum  mechanically,  not  classically. 

( 

Exact  solutions  of  the  Schroedinger  equation  are  possible 
for  some  systems  of  atomic  levels  and  polarizations.  It  is 
particularly  interesting  that  the  stationary  states  of  the  total 
Hamiltonian  of  the  dressed  atom  no  longer  correspond  to  a  fixed, 
well  defined  number  of  photons  in  the  field.  The  number  states 
in>  for  the  pure  radiation  field  are  no  longer  eigenstates  of  the 
total  Hamiltonian.  In  addition,  it  can  be  shown  tnat  stationary 
solutions  exist  which  are  a  superposition  of  one  or  more  of  the 
atomic  states  calculated  in  the  absence  of  fields  and  one  or  more 
of  the  number  states  of  the  field.  This  theory  was  used  to 
predict  the  sideband  intensities  of  well-known  atomic  transitions 
which  were  modulated  by  the  application  of  optical  or  ri  fields 
states.  The  "dressed  state"  theory  predicted  sideband  intensi- 


ties  as  functions  of  the  interaction  energy  (or  field  strength) 


I  very  successfully.  It  is  reasonable  to  inquire  whether  the  same 

•  formalism  could  be  adapted  to  the  analogous  case  in  nuclei. 

:  Maxwell’s  equations,  which  govern  the  interaction  ana  propa- 

j  gation  of  electromagnetic  fields,  are  linear  and  may  be  solved  in 

‘  terms  of  a  Fourier  decomposition.  Magnetization  is  a  more  com- 

r 

f  plex  phenomenon.  and  it  is  expected  that  the  solutions  to  tne 

j  Landau'Li f shi tz  equation  for  obtaining  M(t)  will  contain  a 

!  number  of  harmonics  of  the  frequency  of  the  applied  magnetic 

field.  Each  of  these  harmonics  may  then  be  considerea  to  be  an 
independent  "dressing  field"  in  the  sense  that  Cohen -Tannoud j i 
uses  the  term.  If  done  in  this  way.  the  results  should  be  summed 
over  the  harmonics  which  appear. 

With  this  approach  in  mind.  an  extensive  study  was  mace  of 
the  publications  concerning  the  dressed  atom.  The  most  compre¬ 
hensive  treatise  was  written  by  S.  Haroche  ana  published  in 
Annales  de  Physique  in  1971  in  French.  The  results  of  the  lit¬ 
erature  search  were  condensed  into  an  internal  document  titled 
"The  Interaction  of  the  Dressed  Atom",  copies  ot  which  are  avail¬ 
able  upon  request. 

As  before.  a  better  understanding  of  the  maqne t i za t i on  is 
required  before  this  approach  could  be  adapted.  The  next  sections 
describe  models  of  how  the  direction  of  the  magnetization  varies 
as  a  function  of  time. 

3 )  Development  of  the  Coherent  Switc.h  .Model  of  the  .Magnetization 
Mftl  by  Integrating  the  Landau-Li  f  shi  tz  Ea  ua  t  i  g^n 

As  men tioned  in  the  previous  sections,  both  the  extension  of 

the  early  De.^^aola  and  Collins  model  and  the  dressed  state  tneorv 


of  Cohen-Tannoud j i  are  critically  dependent  upon  the  use  of  a 
good  approximation  to  the  complex  precession  of  the  magnetization 
vector,  Mft)  under  the  influence  of  an  applied  rf  field.  This 
section  describes  an  investigation  of  the  solutions  to  the  Lan¬ 
dau-  Lifshitz  equation  appropriate  to  a  tnin  ferromagnetic  foil 
immersed  in  an  rf  magnetic  field. 

There  are  many  solutions  in  the  literature  to  the  Landau- 

Lifshitz  equation  that  describe  static  domain  walls  for  various 

3  ■’  I  12 

configurations  of  the  magnetization  The  domain  walls  in 

57 

an  Fe  foil  of  the  type  generally  used  in  experiments  are  expec¬ 
ted  to  be  ‘'180°''  walls,  that  is  walls  in  which  the  magnetization 

rotates  through  180°  as  the  wail  is  traversec  from  one  side  to 

3  1 

the  other.  There  are  fewer  models  of  moving  domain  walls  ' 

than  static  walls.  The  precursive  model  of  the  demagnetizing 

field  developed  by  DePaola  and  Collins  to  describe  the  magnetic 

12 

tields  interacting  with  the  nuclear  moment  was  taken  from  Chen 


He  had  described  the  concept  of  the  magnetization  rotating  under 
the  influence  of  the  demagnetizing  field  which  in  turn  is  due  to 
the  magnetic  poles  developed  when  the  applied  magnetic  field 
causes  the  magnetization  to  process  out  of  its  initial  plane. 
Thus  a  moving  domain  wall  is  formed.  This  was  the  mocel  presen- 
ted  by  DePaola  and  Collins*".  On  the  other  hand.  there  is  no 
reason  to  expect  the  foil  to  be  a  single  domain  with  only  one 
moving  domain  wail. 

Experimentally,  it  has  been  found  that  the  static  spectrum 
57 

of  the  Fe  foil  with  no  applied  rf  field  demonstrates  the  usuax 
six  hyperfine  lines.  The  relative  intensities  of  the  parent 


lines  in  this  case  suggest  tnat  the  magne::;  :ield  causing  :ne 
splitting  Ithe  magnetization)  is  con:inea  to  cne  p.ane  cerper.ci- 
cular  to  the  transmission  axis.  that  is.  t.ne  plane  o:  t.ne  :oi.. 
When  an  rl  magnetic  tielh  vs  applied.  tne  spectrum  cecomes  muc.n 
more  complex.  There  are  the  usual  six  hvperrine  lines  :pa.-en: 
lines)  plus  six  systems  of  sidenands.  tnat  is.  apsorpticn  *ines 
displaced  from  the  parent  line  bv  trequencies  cor respcnc  in g  to  an 
integer  (positive  or  negative;  numner  times  tne  3pp.i.ed  :.'e- 
quency.  The  relative  intensities  cl  cne  parent  lines  continue  to 
suggest  that  the  magnetic  lield  remains  contined  to  the  plane 
perpendicular  to  the  transmission  axis. 

Thus,  the  internal  magnetic  structure  o:  tne  '  Fe  :oi.  is 
considered  to  be  a  composite  of  many  p  i  a t e . e t - s na pea  mitro- 


crystals  with  the  platelets  aligned  more  or  less  paralle. 
Che  plane  of  the  foil.  Demagnet  i  c ng  tie.os  cevelop  in 
volume  of  cne  platelets  due  to  tne  appearance  c:  maine’.t  : 


on  the  surfaces.  Each  microcrvstal  or  plate. e*:  is  tor.siaerec  tc 
be  a  single  domain.  In  tne  conerent  swittnin;  mcce..  t ?  ma;ne- 
t i za  1 1  on  at  ail  points  wi t  n i n  tne  ccma . n  '  t  ■  a  ■  e  s  in  p  na  s  e .  ■ n  a • 


together  as  a  unit  witncut  'ne 


t  m  a  1  .n  w  a  .  .  w  r.  e  ^  o 


rotation  is  concentrated  witnin  tne  wa  .  . 


are  not  identica..  tnere  are  sister,  tin  -  i?i 


latnatiti^:  :n 


between  the  platelets,  put  insice  tne  p.a'e.e^,  "ne  ma  sn*^  "  ;  c  a  '  ;  :  r. 
is  uniform.  se-cn o  m.ode  wnic.n  in-..ces  "no  ^::o:t;  ::  tn^ 

traversal  ot  a  domain  wa . .  tnrt  .gn  a  n.c.ear  ;:"e  w.  .  .  c-  ;;n 
sidered  in  section  4. 


Equation  (D;  aestripec  ‘  n  m :  n 


ence  of  tne  four  magne': 


a  .  an:  : 


previously.  The  problem  was  solved  using  a  combination  of  spher- 
ical  polar  coordinates  and  the  Cartesian  components  of  Ml t )  as 
shown  in  Figure  8.  The  Cartesian  components  of  Mlt)  are 


M  =  M  sin  8  cos  <t> 

X  s 

M  =M  sin^sind 

y  s 

M  =  M  cos 
z  s 

These  four  fields,  the  applied  field, 
tield,  the  anisotropy  field,  and  the  exchange 
considered  as  part  of  this  model. 


C  3  4  a  J 
i34b) 

I  34c  1 

the  magnetostatic 
field,  will  each  be 


H  l  t 
app 


!1 .  sin  (j ,  t 

X 


A)  The  Applied  Field.  Consistent  with  the  initial  condi¬ 
tions.  the  applied  magnetic  field  is  written  as 

135) 

3)  The  Magnetostatic  Field.  The  .Magnetostatic  field  is 
given  by 

B  =  li^lH  ^  M)  136) 

wnere  H  is  given  by  the  solution  to  the  Poisson  equation 


H  =  1 1  /  4-!r  )  grad 


f  di_v  M  .  r  1:1 

^ii  -'s  "ii 


137) 


IS  denoted  bv  r...  The  first  intedrai  is  taken 
1  1 


The  distance  between  the  point  of  integration,  i  surrounded  by 
the  volume  of  integration  dV .  and  the  point  j.  at  which  .4  is  to 
he  evaluated. 

throughout  the  entire  volume  of  the  magnetized  bodv.  The  second 
integral  is  taicen  over  the  entire  suriace  of  the  magnetized  body, 
n  being  tne  inwardly  directed  unit  vector  normal  to  the  suriace. 

The  magnetostatic  field  is  best  considered  on  a  case*  by- 
case  oasis.  The  divergence  of  a  vector  field  is  a  measure  of  the 
'sources  or  "sinks"  of  that  vector  within  a  volume.  Since  tne 


magnitude  of  M  remains  the  same  througnout  the  boay.  div  M  =  0. 
The  first  integral  resembles  the  volume  integration  of  a  charge 
density  equal  to  div  M  which  vanishes  throughout  the  volume. 


The  second  integral  resembles  an  integration  of  a  surface 
charge  density  over  a  surface.  The  non-zero  magnetic  pole 
density  arises  from  the  discontinuity  of  M  at  the  surface  and  the 
continuity  boundary  requirement  on  B  and  H.  If  M  has  a  component 
normal  to  the  surface,  the  magnetic  pole  density  is  CM*nj  on  the 
surface.  As  the  applied  magnetic  field  increases  parallel  to 
the  X  axis,  the  magnetization  begins  to  precess  about  the  applied 
field,  and  MCtl  moves  out  of  the  xy  plane.  A  magnetic  pole 
density  develops  on  the  faces  of  the  foil  parallel  to  the  xy 
plane  (boundary  conditions  on  3  and  H).  These  magnetic  poles, 
whose  density  is  given  by  *  ( t 1 ,  are  the  source  of  the  magneto- 
static  field,  which  then  simplifies  to  the  single  term, 

H,Ct)  =  -  Ml  =  -  M  cos  e  .  (38) 
d  s  2 

Cl  The  Anisotropy  Field.  The  ferromagnetic  foil  is  assumed 
to  consist  of  many  small  microcrystals  shaped  like  platelets  or 
"microfoils".  This  is  consistent  with  t.he  amplitude  d  i  s  t  r  i  pu  t  ion 
of  the  observed  six-line  spectrum  in  the  absence  of  any  applied 
field  which  is  a  rather  sensitive  function  of  the  average 
orientation  of  the  magnetization  of  each  domain  in  the  sample. 

A  magnetic  material  is  said  to  have  an  easy  or  preferred 
direction  of  magnetization  if  a  minimum  energy  state  is  ootained 
by  having  the  magnetization  lying  along  this  easy  direction.  In 
a  uniaxial  material  (the  approximation  used  for  each  platelet  or 
the  ferromagnetic  foil),  any  deviation  of  (  t  j  away  from  the  easy 


direction  by  an  angle  results  in  an  increase  in  the  energy 

density  which  may  be  accurately  represented^  by  the  expression 
K  sin‘^(<^-0  ].  The  anisotropy  constant  K  has  units  of  energy/vo- 
lume . 

3 

The  effect  of  anisotropy  is  often  represented  by  an 

equivalent  Cbut  fictitious)  magnetic  field  which  lies  along  the 

easy  direction  of  the  crystal.  The  exchange  field  and  the 

anisotropy  field  are  derived  in  the  same  way;  the  expression  for 

the  inte. .  ac  t  ion  energy  is  derived  from  basic  principles  and  then 

a  fictitious  magnetic  field  is  deduced  which  yields  the  same 

functional  expression  for  the  energy,  based  on  the  relation, 

E.  ,  =  -Cl/2)  u  M-H.  139) 

1  n  t  o 

A  unit  vector  in  the  direction  of  the  easy  axis,  as  shown  in 
Figure  8,  is  designated  Thus  the  anisotropy  field  is  given 

by 

=  Q  CM*e^,)  e^,  =  Q  sin  8  e., ,  C40) 

where 

A  /^ 

e  ,  =  cos  d  e  -*■  sin  0  e  141) 

X  ox  o  y 

2 

and  where  Q  =  C 2K  /u  M  !  is  a  parameter  of  the  crystal.  The 

u  o  s 

angle  is  the  acute  angle  between  the  easy  axis  anc  the  ;c  axis 
of  the  foil  coordinates. 


D)  The  Exchange  Field.  Since  t.ne  magnet  icat  ion  at  every 
point  within  the  domain  precesses  in  the  sa.me  manner.  that  is, 
in  the  domain  where 


MCr.t)  =  MCr+dr.t) 

the  exchange  field  is  zero.  .However,  in  materials 
depends  on  the  position  of  the  point  of  oDse r v a 1 1 on . 


wne  r  e 


such  as 


in  a 


28 


Landau'Li  f  shi  tz  wall,  the  "exchange  field"  !  i  c  1 1 1 1  ous  J  :s  given 
by 

Hir.tj  =  i2A/p  M“J  Mfr.t)  '4J; 

o  s 

where  A  is  a  constant  of  the  material. 

The  total  magnetic  field  F  in  equation  i2j  is  the  sum  ot  the 
fields  given  in  equations  (35).  (38)  and  (40;: 


(H 


app 


u 

■*d 


H  ) 

a 


44  , 


The 

reduced 
terms  of 
sin 


equations  of  motion  of  the  components  of  M',  t  ;  can  be 
and  combined  to  yield  equations  for  dii/dt  and  d^/dt  in 
the  components  of  F. 

8  (dd/dtj  =  -  I  v  I F  [  (.\/M  )sin  i  ♦  cos  t  cos  6\ 

X  s 

I>IF  [  (A/.M  )cos  *  -  sin  *  cos  8]  (45) 

y  s 

I  y  !  sin  8 


(  d8  /  d  t  )  =  -  !  >  '  ‘  2  in  4)  *  (  A/M^  )  cos  l>  cos  8] 

•^!>IF  iCosii*iA/M)sinif>cos8'’  (46) 

y  s 

-  !  y  !  F  ^  (  A  M  ^  ;  sin  8 

The  components  ot  the  total  magnetic  lieid  F  are  given  below. 


c  = 

X 

-  u  H, 
o  1 

sin  w.t  *  u  M  Q 

1  0  3 

sin  8  c  o  s '  p  -  t  _  j  cos 

Q 

,  4  ~  a  : 

c  = 

‘  y 

♦  11  .M 
o  s 

Sin  8  cos  [  t  -  i ^  } 

S  ITi 

■s 

,  4':’b  ; 

P  r 

*T 

-  U  .M 

o  s 

cos  8 

1  4  ■('  c  ' 

These 

are  f 

i r  s  t  order  differ 

en  1 1  a  1  s  in  t  i  .me  out 

V  e  r  ■/ 

comp . 1 

cated  functions  ot  the  angles  8  and  i.  .A  crcira.m  was  written  in 
FORTRAN  for  the  .Hewlett-Packard  Series  3000 '560  'I'tmputer  in  wnicn 
the  angles  8  and  <6  were  Duilt  up  in  small  steps  tro.m  tne  initial 
conditions  to  obtain  the  magnetization  .Mit:.  An  enamo.e  solution 
is  shown  in  Figure  SA . 


The  motion  ot  Ml t  '  seen  in  Figure  BA  is  gui 
with  the  expectations  reached  by  purely  qua  i  1 1  a  ■:  i e 
t  =  0.  tne  magnetization  lies  in  tne  piane  o:  tr.e  :o 
radians  along  the  easy  axis  at  an  initial  ang^e  ^ 
be  0.35  radian  in  this  example.  The  initial  angles 
solution  parameters.  As  the  magnetostatic  iield  me 
tne  precession  of  M; t }  out  ot  tne  xy  piane,  tne  e 
magnetostatic  rield  I  in  the  z  direction  com?  1  e :  e  * 
motion  of  MitJ,  causing  it  to  rotate  in  t.ne  xv 
until  It  nears  tne  easv  axis.  Alter  sin  uj.t  reacn 
and  oegins  to  decrease,  tne  .M^  begins  to  decrease  a 
to  the  xy  piane.  The  magnetization  now  .les  a .on; 
Put  anti  parallel  to  its  initial  direction.  'Tne  m 
proceeds  in  a  similar  manner  tnrougn  tne  sec  one  na. 
applied  field  (aitnougn  tne  motion  is  rot  a  tec  aoc 
through  130"^  from  its  initial  position.. 

The  angle  i  t  )  as  a  function  of  time  is  also  s 
9.A .  The  oscillations  of  .M'  t  ;  aoout  -ne  eoui.io 
between  tne  easy  direction  and  tne  direction  c:  t ? 
are  clear  iv  s.nown.  Wnen  M  c  ’  switc.oes.  it  is  :a 
equiliorium  position.  Tne  second  "erm  ecua":.,". 
a  r  e  s  t  o  r  1  .n  d  force  w.n  i  c  .n  c  r  i  n  c  o  M  .tic  re  n  e  a  r  .  e  r  a 
equilibrium  cireocicn.  The  equi.icrium  d.recticr. 
oy  the  local  minimum  o:  ■:ne  sum  c:  'ne  in-eracti: 
tne  magnetization  witn  n  e  an  i  so  t  r  c  p  v  rie.t  an  c 
field.  I:  tne  applied  fie.o  were  verv  cm:,  r.  :. 

r  e  s  t  o  r  1  n  d  force  would  c  r  i  n  c  t m  a  d  r.  ?  t  ;  c  a  m  ■  n  car 


cons  1 s  tent 
umen  t  s .  At 

t  9  ==  r'  2 

assumed  tc 


menaces  tne 
ne  rapid].'/ 
Its  maxi  mum 


applied  field 


1  j  u  r  e  9  b  .  w  r. .  c  n 


CO  confirm  this  initially  surprising  result. 

Figures  lOA  and  lOB  show  the  dependence  ot  tne  motion  oi 
Ml  t ']  on  the  initial  value  oi  the  angle  S  for  the  values  of  d 

o 

equal  to  91°  ana  89°,  respectively.  From  these  ligures  it  is 

evident  that  the  starting  value  of  6  determines  the  direction  oi 

o 

the  initial  rotation  only.  After  tne  initial  eye  1 e .  the  Peha vior 

of  Ml'tl  is  the  same  in  all  three  cases.  Since  the  magnetization 

is  expected  to  lie  in  the  plane  of  the  foil  fwhich  minimizes  the 

magnetostatic  energy]  only  values  near  90  can  be  realistic.  The 

maximum  expected  deviation  from  90°  would  be  on  the  order  of 

6d  (thickness  of  the  f oi  1 1  inear  dimension  of  the  foil) 

o 

~  5  X  10  ^  m  /  1  X  10  “  m  ~  0.03  degree 

For  this  reason.  an  initial  value  of  d  ''  90°  is  considered 

o 

r easonabl e . 

Finally,  parameter  studies  have  snown  tnat  for  small  values 

of  .H,  and/or  certain  values  of  b  ,  the  direction  of  magnetization 
1  o 

in  the  foil  does  not  switch.  This  is  consistent  with  the  picture 
above  in  chat  the  applied  field  never  develops  enoug.*!  interaction 
energy  to  overcome  the  effective  anisotropy  energy  barrier,  ^  . 

The  most  important  output  of  this  model  is  the  decompos i t i on 
of  sin  .^it;  and  cos  (  t  j  in  terms  of  harmonics  of  the  applied 
field  frequency.  The  amplitudes  of  the  harmonics  then  beco.me  tne 
input  to  tne  program  tor  solutions  to  the  time  dependent 
Scnroedinger  equation  from  which  sideoand  intensities  are 
computed  according  to  one  or  another  of  ‘he  approaches  oeing 
discussed.  The  amplitudes  of  the  harmonic  decomocs i t i ons  of  cos 
t)[  t  i  and  sin  ii  i  t  ;  versus  the  order  ot  the  harmonic  are  snown  in 


Figure  11  for  a  typical  sec  of  input  parameters. 

Since  tPis  type  of  numerical  integration  of  the  Lancau- 
Lifshitz  equation  must  be  conds ider ed  "exact"  to  witnin  the 
limits  of  numerical  methods.  its  principal  shortcoming  is  tnat 
coherence  is  limited  to  a  single  domain.  The  possible  importance 
of  domain  walls  is  considered  next.  That  is.  this  picture 
requires  the  magnetization  of  the  entire  domain  to  switcn  without 
the  presence  of  a  moving  domain  wail. 


In  Che  coherent  switching  model,  magnetization  within  the 
domain  is  required  to  rotate  uniformly  (equation  42).  It  is  also 
possible  fo"  the  direction  of  the  magnetization  to  change  by 
means  of  passage  of  a  domain  wall  through  the  material. 

There  are  several  different  kinds  of  domain  walls.  for 
example  the  90*^  wall  where  the  direction  changes  by  90°.  and  the 
180°  wall  between  two  domains  wnose  magnetizations  are 
an t i pa r a  1 1  el .  An  example  containing  a  180°  wall  is  snown  in 
Figure  12. 

The  position  and  shape  of  a  domain  wall  is  determined  py  tne 
minimization  of  all  the  interaction  energies.  In  tne  absence  o : 
an  applied  magnetic  field,  tne  domain  walls  remain  static.  Vhen 
a  magnetic  field  is  applied.  the  domai.n  wail  moves  in  sued  a  wav 
that  the  volume  of  the  domain  wnose  magnetization  is  parallel  :c 
the  applied  field  increases.  The  static  and  moving  domain  walls 
will  be  discussed  next. 

A)  The  Static  Domain  Wail. 

Within  the  domain  wall,  the  direction  of  M  rotates  wnile  its 
remains  fixed  at  the  saturation  value  M  . 


raagn i tude 


The  fundamental  equation  of  microraagnet ics  states  that  the 
torque.  given  by  the  vector  product  of  the  magnetization  ano  the 
total  magnetic  field,  must  vanish  at  every  point  in  a  medium 
which  is  in  magnetostatic  equilibrium.  That  is, 

(M  X  F)  =  0.  148) 

For  an  infinitely  long  sample.  the  H  fielc  calcauiated  from 
equation  (37]  is  neglible.  Since  the  magnitude  of  M  remains 
constant,  the  divergence  vanishes.  There  is  no  normal  component 
of  M  at  the  surface  of  the  sample,  except  either  on  faces  which 
are  far  removed  in  the  z  direction  from  the  point  of  observation 
or  at  the  very  small  area  where  the  comain  wall  meets  the 
surface.  Therefore,  the  contribution  to  F  from  the  magnetostatic 
field  is  simply  by  equation  (37). 

However,  since  M  is  not  uniform  within  the  wall,  there  is  a 
contribution  from  the  exchange  field.  The  easy  direction  of 
magnetization  is  taken  to  be  parallel  to  the  z  axis.  The 
components  of  F  are 

F  =  -C2K  /M")  M  +  (2A/M^)  +  u  M  (49a] 

X  USX  s  xox 

F  =  (2A/M^]  v^M  ^  w  M  .  r49b) 

2  S  Z  O  2 

In  this  case  there  is  assumeq  to  be  no  appliec  magnetic 

field.  SO  that  H  =  0.  The  equation  determining  the  spatial 

app 

dependence  of  M  is  obtained  by  substituting  equations  [49a]  and 
(49b)  into  equation  (48), 


which  is  an  ordinary  differential  equation  since  M  is  a  function 

of  y  only.  The  components  of  M  may  be  written 

M=Msin^  C5la') 

X  s 

M  =  M  cos  8  ,  1 5 lb) 

z  s 

and  the  solution  is  found  to  be 

±  tan(5/2)  =  exp  Ly/A)  [52) 

where  A  is  the  well  known  wall-width  parameter, 

A  =  (A/K  .  [53) 

u 


One  should  note  that  the  angle  8  is  that  angle  8  defined 
Figure  12. 

The  choice  of  sign  in  equation  [52)  indicates  that 
domain  wall  may  have  either  a  clockwise  or  an  an t i c 1 ocKwi se  sc 
sense  associated  with  the  rotation  of  M  in  the  xz  plane, 
components  of  M  are 


in 


t  he 
r  ew- 
Thi 


M  =  -  M  tanh  [y/A)  [54aj 

X  s 

M  =  +  M  sech  [y/A)  .  [54b) 

z  s 


3)  The  Moving  Domain  Wall. 

In  1935  Landau  and  Lifshitz  applied  their  equation  [equation 
[2)]  to  the  case  of  the  motion  of  the  static  domain  wail 
structure  described  above.  Thev  assumed  that,  although  the  wall 
was  moving  due  to  the  influence  o:  the  applied  magnetic  field, 
the  wall  structure  remained  almost  identical  to  its  previous 
static  structure  [equations  [54a)  and  L54bj].  That  is.  they 
assumed  a  "rigid  wall  structure". 

If  the  wall  moves  as  a  rigid  structure.  then  the  spatial 
derivatives  of  M  are  defined.  If  the  wall  moves  at  a  constant 


velocity,  the  time  derivatives  are  defined  as  well.  The  problem 
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becomes  algebraic  and  may  be  solved  ac  anv  convenient  point,  such 

as  the  center  of  the  domain  wall.  At  the  center  oi  the  aomain 

wall,  M  =  M  and  M  =0.  The  other  assumption  wnich  loilows, 
X  s  z 

once  a  rigid  wall  model  is  adopted,  is  that  the  anisotropy  field 
and  the  exchange  field  remain  antiparallel  to  M  in  the  moving 
wall  as  in  the  stationary  wall.  Then  only  the  applied  field  must 
be  considered. 

The  domain  on  the  left  in  Figure  13A  is  magnetized  in  the 
same  direction  as  the  applied  field  and  consequently  grows 
through  the  motion  of  the  wall  to  the  rignt  with  velocity  v^. 
This  figure  only  shows  M  in  the  two  domains  and  at  the  center  of 
the  wall,  assuming  a  positive  screw-sense  of  the  rotation.  Here 
only  B  enters  into  equation  (.37): 

“7 

CdM  /dt )  =  ! V  IM  B  1 55a ) 

y  s  z 

(dM,/dt)  =  ‘vIAB  .  ''55b' 

z  z 

The  damping  constant,  A.  can  only  be  cetermined 
experimentally  by  ferromagnetic  resonance  on  single 

3 

crystal,  spherically-shaped  samples.  In  general"",  the  di.mension- 
less  ratio  (A/. Ms)  has  been  found  to  he  on  the  order  ol  I'j  ^  or 
less.  Since  (W/M^)  "  10  equation  !'55a)  is  the  more  important 
of  the  two.  This  equation  indicates  that  tne  proposal  that  the 
wall  can  move  forward  under  the  influence  of  the  apclieo  tielc 
must  involve  tne  vector  M  developing  a  component  in  tne  airection 
of  motion.  It  is  not  possible  to  assume  tnat  tne  wall  is  reallv 
rigid  and  maintains  exactly  the  same  form  which  it  has  wnen 
stationary.  As  shown  in  Figure  13B.  M  must  tilt  bv  the  angle  d 
to  satisfy  equation  (55a). 

At  the  center  of  the  wall,  t.ne  components  of  .M  are 


C56a  J 


M  =  M  cos  <t> 

X  s 

M=Msin0  (5Gbj 

VS 

=  0  .  i56c) 


In 

the 

original 

development 

of  the  model,  Lanaau 

and 

Li f  shi tz 

continued  to 

assume  that 

the  contributions  to  F 

due 

to 

the  exchange 

and  anisotropy  fields 

remain  antiparallel  to 

M  _ 

as 

when  0=0.  The  component  is  assumed  to  have  the  same  torm  as 


gi7en  in  equation  C54b).  Then  equation  (37)  is  sol7ed  using  only 
the  magnetostatic  field  and  the  applied  field. 


-(dM^/dt'J  = 

l>l 

sin  0 

(B  M 
z  s 

■  A/j  M 
o  s 

cos 

0  s in  0 ]  C  57a  ) 

■"CdM  /dtj  = 
y 

1  y  1 

cos  0 

[B  M 
z  s 

-  Ay  M 

0  s 

cos 

0  s  1  n  0  ]  (  5  7  D  ) 

*(dM  /dt)  = 
z 

1  y  1 

[pnl 

0  S 

cos  0 

sin  0  + 

AB  : 

1  (57c) 

From  this,  one 

sees 

Chat 

domain 

wall  mot 

ion 

may  be  visualized  as 

follows.-  as  the  wall  mo7es  forward.  the  7ector  M  rotates  aoout 

y.  Since  M  has  constant  magnitude  and  *  is  constant  for  a 

constant  7elocity,  equa  tions  C57a)  and  (57bj  are  zero.  The  angle 

0  is  related  to  the  applied  field  through 

s  i  n  2  0  =  !,  2  3  /  u  A )  .  !'  5  3 

z  o 

Since  a  rigid  wall  nas  been  assumec.  at  the  center  o:  the 

wa  1 1 

cM  /dt  =  !  dM  /3y)(dy/dtJ  =  I  .M  /d)  v  ,  '53: 

z  z  s  y 

wnich  wnen  combined  witn  (57c)  leads  to  an  expression  :or  v  .. 

2 

7  =  (!7!d/M)[AB  sin0cos0l  .  -  ~z'j 

V  S  Z  o  3 

This  gi7es  the  picture  of  a  domain  wail  mo7ing  torwa.c  at 

constant  7elocitv  7  under  the  influence  of  a  constant  acc.iec 

y 

field  S_ .  Combining  equations  (58)  and  (60).  t.ne  velocity  ic 


given  by 


=  ClylAM^/Ajfl  A“/mJ] 


(SI  ) 


In  conducting  materials  where  eddy  currents  complicate  the 
calculations,  it  is  possible  to  derive  an  expression  for  the 
domain  wall  velocity  in  terms  of  the  product  of  a  domain  wall 


mobility  and  the  applied  field,  that  is. 

7  =  u  B 

y  w  z 

The  domain  wall  mobility  is  given  by 

h]  ^ 

w  os 


[62  j 


[631 


where  <r  is  the  conductivity  and  h  is  the  thickness  of  the 
material. 

From  this  domain  wall  velocity,  a  tilt  angle  a  coulc  be 
deduced.  If  <t>  developed  to  approximately  45°,  its  maxi.mum  pos* 
sible  value.  a  significant  amount  of  the  magnetization  field 
might  be  transmitted  to  an  adjacent  material.  .However.  it  was 
found  t.hat  t.he  passage  of  a  wall  at  the  surface  transferrec  only 


on  the  order  of  0.1  to  10.0  H 


wnich  was  calculated  to  oe  at 


tne  most  20  gauss.  Excitation  of  single  non-magnetic  srainless 

57 

steel  foils  enriched  with  Fe  at  tnis  levei  failec  to  snow 
sidebands.  Thus  it  was  concluaed  tnat  tne  wall  passage  a* one 
could  not  add  significantly  to  the  transfer  o:  an  effect  from  a 


magnetic  to  a  non-magnetic  foil  at  tne  level  of  magnituce  seen  in 


exper iment  s . 

In  general,  computer  analyses  of  the  static  oomain  wall 
equations  in  conductors  yield  only  solutions  in  whicn  tne  magne¬ 
tization  remains  parallel  to  the  surface.  .A.  .Anaroni  nas  mace 
extensive  analyses  of  domain  walls  in  conauctint  materials.  and 
his  solution  for  a  uniaxial  material  is  snown  in  ."igure  14. 


37 


These  analyses  reinforce  the  conclusion  that  the  effects  of 
moving  domain  walls  in  conductors  are  insufficient  to  transfer 
excitation  to  the  non-magnetic  foil.  The  problem  of  the  transfer 
of  sidebands  from  magnetic  to  non-magnetic  layers  lies  at  the 
focus  of  evolving  concepts  for  the  next  year,  but  it  must  be 
recognized  that  such  transfer  leads  to  a  relatively  small  effect, 
anyway.  The  largest  effect,  by  an  order -of -magni tude ,  is  the 
development  of  sidebands  in  a  single  magnetic  foil  and  this  seems 
to  be  able  to  be  described  by  the  model  synthesized  from  the 
previous  approaches  as  described  in  the  following  section. 

5J  Solutions  to  the  Time-dependent  _Schrogdipqer  Sguacioh  ^ 

Higher  Level  of  Approximation 

A  more  tractable  approach  to  the  solution  of  the  time  depen¬ 
dent  Schroedinger  equation  discussed  earlier  was  suggested  by 
Professor  Paul  Berman  of  New  York  University.  It  is  similar  to 
some  developed  earlier  but  proves  useful  in  more  general  cases 
for  single  foils. 

In  con  t  inued  collaborati  on  with  Professor  Berman ,  t  .h  i  s 
approach  has  been  programmed  for  the  Hewlett-Packard  Series 
9000/560  Computer  and  is  in  the  process  of  checkout.  This 
approach  combines  the  exact  solution  of  the  time  dependent 
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Schroedinger  equation  for  the  system  of  the  Fe  nucleus  in  a 
ferromagnetic  foil  with  perturbation  theory  solution  of  the 

absorption  of  the  gamma-ray  photon. 

Figure  7  shows  the  energy  levels  and  angular  momentum 

57 

assignments  of  the  Mossbauer  hyperfine  transitions  in  Fe .  Be¬ 
fore  detailing  the  calculations.  the  matrix  indices  of  the  va- 


nous  states  ct  the  nucleus  must  Jetinea.  To  avoio  the  un¬ 
necessary  complication  ot  carrying  octh  J  anc  M,  is  subscripts 
denoting  a  specific  level.  eacn  magnet;;  suoieve.  tne  jreunc 


state  and  excited  state  has  seen  assignee 


a  single  i n  e  e  x  a , 


f  ol lows : 
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The  total  Hamiltonian  0:  tne  re  nucleus  ;  r.  a  immerse; 

in  an  rf  magnetic  lielc  is  given  pv 

H  =  *  V, I : .  ♦  V,  t  64 

01 

where  H  is  tne  nuclear  Hamiltonian  wi  t.n  V  :■  =  D  anc  no  caama 

o 

ray  interaction.  The  matrix  H.  is 

r  0  0  '  j  1'  0 

I  0  'j  !  'J  : 


In  Eg.  1641.  V  ,  t  i  is  the  octer.tia*  tne  1  r.  t  e  r  a  ;  t  1  ;  .n 
magnetization  M',  t  with  -.ne  nuc.eus  navin;  magnet.:  miTen"  - 

\I  I  r  I  r  -  , ,  y  1  ^  I  :  ;  -  _  ,  M  .  ■ 


V^i  t  J 


where  g^  is  the  gyrcmagneti:  ratio  ;•  *  the  nuc.e 
is  the  nuclear  magneton.  and  J  is  the  7e'':o 
momentum  operator. 


-  a  .  a  h  ;  u  .  a  r 


1^*  •  ^ ^  *"%  •*«•*«•* 


.N  _ 
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The  matrix  V^(t)  is  given  by 


V^C  t) 


• 

q  M 
g  z 

3/ 

-i 

1 

0 

0 

0 

g  M 

9/ 

zj  ^ 

0 

0 

0 

0 

0 

:  3g  M 

1  '“e  z 

1 

3g^M. 

0 

0 

0 

0 

i 

Sg^M. 

0 

0 

0 

1  0 

1 

■®e”z 

3ggM. 

0 

0 

i  0 

t 

0 

3geM, 

-3g^M 

is  the 

same 

as  shown 

in  Figure 

8. 

If  Mlt) 

C67) 


component,  the  diagonal  elements  of  V^(t)  vanish.  The  compo¬ 


nents  of  Mltj  in  the  plane  of  the  foil  are 

U  -  M  1*1  i  M  f*-1  -  U 


=  M  itl  +  i  M  (t)  =  M  C68) 

*  X  y  s 

The  time  dependence  ol  the  angle  *(t)  is  determined  from  the 
techniques  described  in  Sections  3  and  4. 

The  coupling  of  the  gamma  ray  field  to  the  nucleus  is  given 
bv  The  selection  rules  for  nuclear  magnetic  dipole  trans¬ 


itions  (Ml  transitions)  are  A J  =  0 ,  1  and 


tl.O.  The 


matrix  VTl'tj  has  the  torm 


V  !  t  •  = 


where  K 


0 

1 
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t  he 

matrix  element  ol 

t  he 

gamma 

C69) 


operator . 


■'ilV  cos  w  tlk/ 
>  y 


The  time  lepenaent  Schroedinger  equation 
i  h  (' da  /  d  t  )  =  H  a 
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is  written  for  this  system  as 

iM  (da/dt)  =  +  VAt)  +  V„it))  a  . 

o  1  2 


(72) 


The  problem  can  be  divided  into  two  parts.  First,  the  time 

dependent  Schroedinger  equation  is  solved  exactly  for  the  Hamil" 

tonian,  H  =  H  V,Ct).  Then  the  solutions  are  used  to  treat 
o  1 

V2Ct)  as  a  perturbation  of  the  system. 

First,  the  state  of  the  system  is  transformed  into  the 


interaction  representation, 

r  .  -ant/Yi)  - 

a  -  e  o  D 


C73) 


Explicitly, 

for  k  =  1,2  C ground  state) 

for  k  =  3, 4,5,6  [excited  state). 
When  equation  C73)  is  substituted  into  equation  (72),  the 
equation  for  (db/dt)  becomes 

ih  (db/dt)  =  (V’  ^  V;)  b  (74) 


where  the  V  matrix  is  defined  by 

V-  -  oCiH^t/h)  ^-(iH„t/M) 
V  -  e  o  V  e  o 


or 


V’  .  .  =  V.  .  e''  ^“i  i  ^ 

13  1 3 

‘**i  3  "  ^  “i  '  E  j  )  /W 


(75) 

(7Sa) 

(76b) 


Since  V^(t)  does  not  have  matrix  elements  between  the  ground 
state  and  the  excited  state,  V’^(t)  =  V^(t). 

However,  V2(t)  does  have  matrix  elements  linking  the  ground 
state  to  the  excited  state.  The  non-zero  matrix  elements  of  VI 


are 


■  "13  =  '''2>*31 


(77a  ) 


(  V  ’  =  K  e'  ^*^0^ 

^2^14  ‘^14  ®  ° 


(v;) 


2 '41 


(77b) 


The  gamma  ray  interaction  may  be  written  as 

V.  =  V  cos  «  t  =  Cl/2)  V  +  e‘^“y'^)  .  C78] 

2  >  V  y 

The  rotating  wave  approximation  consists  of  neglecting  all  of  the 


terms  which  contain  the  sum  of 

the  frequencies, 

+  (a;  + 

y 

u>  )  ,  in 
0 

the  exponential  time  dependence. 

In  this  approximation. 

the  non- 

zero  matrix  elements  of 

become 

<'^>13 

^i At  . 
e 
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t  79b ) 
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i  At 
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<''2’S1 

C79c) 

'''2>24  = 

'\)24 

i  At 
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CV’ Igj 

(79e) 

<''^26  = 

<\’2S 

i  At 

e  = 

CVj )g2 

1^91  ) 

where  the  angular 

frequency  L  is 

defined  by 

&  -  U) 

y 

Ck^ 

0 

18O) 

The  matrix  UCt)  is  defined  as  the  V'2  matrix  in  the  rotating 
wave  approximation  as  shown  in  equations  C79a)  through  (79f) 
above.  In  the  rotating  wave  approximation.  equation  ('731  is 


r ewr i t  ten . 


iH  Cdb/dt]  =  CVj  +  U)  b 


(81 ) 


Furthermore,  the  vector  b  may  be  written  as  the  matrix  product 


b  =  exp(( 


r 

i/hll  Vj(t’ 

Jo 


)dt ' ]  c  5  TC  t )  c  . 


(82] 


The  equation  of  motion  for  the  matrix  T(t]  is 


with 


iM  (dT(t]/dt]  =  Vj  T(t)  (83] 

the  initial  condition,  TCt=0]  =  1.  tne  unit  matrix. 

The  equation  of  motion  of  the  vector  c  is  now  given  by 

iM(dc/dt]=U’ft]c  (84] 


where  U’(t]  is  defined  by 

U’(t]  =  T"^  U(t]  T  .  (85] 

The  matrix  (J  ’  ( t  ]  has  the  same  general  form  as  shown  above 
in  equations  C79a]  through  (79f).  The  ground  states  and  the 
excited  states  are  now  mixed  as  a  result  of  both  and  The 

relaxation  effects  of  the  natural  decay  of  the  excited  state  are 
now  added  into  equation  (84], 

(dc/dt]  =  ( -i/h]  U’ ( t  ]  c  -  ( 1/2  ]:^  c  (86) 

where  the  matrix  y,  is  defined 
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The  quantity  y  has  the  dimensions  (sec  and  is  defined  as 
inverse  of  the  natural  lifetime  of  the  excited  state.  y  =  ( 
The  net  absorption  of  the  gamma  rays  is  proportional  to 


the 


1  /  r  }  . 


6 

Z 


a  (t]l  for  times 


t 


>  > 


( 2 / V  J .  However  , 


9  “  9  9 

I  laCt]|“=  I  Ib(c3!"=  I  icltj!“ 

3]S  o  ^  ^  ^ 

m=3  m=3 


(80  ) 


since  and  are  Hermitian.  The  components  o!  c(tj  are  given 
by 


c,Ct3  =  C-i/h)  I 


k  =  : 


r.-, 

1  Jo  ‘ 


For  t  >>  (2/-y),  the  absorption  is  proportional  to 


<  E  ic,  it] r> 
k  =  3 


where  the  time  average  is  taken  over  the  time  response  ot  tne 
detector.  This  completes  the  formal  solution. 

If  the  decay  of  the  excited  state  into  the  ground  state  is 
neglected,  then  to  the  zeroth  order  in  V^it). 

Cjj  ^  ^ 

since  ^2^^^  does  not  couple  the  ground  states.  1  ana  2,  in  tne 

first  order.  Furthermore.  since  H  ana  V,  are  Hermitian,  the 

o  1 

amplitude  Cj^(t)  for  the  ground  states  may  be  written 

CkCt)  =  ^ 

The  ground  state  splitting  is  small  compared  to  tne  thermal 
energy  at  room  temperature,  so  that  the  ground  states  are  equally 
populated.  Thus. 

aj^CO)  =  Cl/v/2)  for  k  =  1.2  02) 


The  absorption  spectrum  then  becomes  the  result  ot  six 
computational  steps. 

Stgp  i:  To  determine  the  V^lt)  matrix  given  by  equations  07). 

equation  (68)  may  be  expanded  to 

M  =  M  (  cos  i*(t)  t  1  sin  41  t  )  }  03) 

s 


The  output  ot  the  program  written  tor  the  coherent 


swi tching 


model  described  above  includes  the  Fourier  decomposition 

amplitudes  A  ,  B  ,  A’,  and  B’.  where 
n  n  n  n 

cos  (*Ct)  =  H  f A  cos  nw^t  +  B  sin  nw.tl  (94) 

n  In  1 

n  =  0 

CO 

sin  li(t)  =  I  fA’  cos  nw,  t  -►  B'  sin  nu.tl  .  (95) 

n  =  0  "  In  1 

These  amplitudes  may  be  used  as  inputs  which  specify  V^(t). 

SleB  2:  Step  2  is  the  numerical  solution  of  equation  (83)  for 

the  matrix  T(t)  subject  to  the  initial  condition,  T(t=0)  =  1. 

Equation  (83)  is  a  first  order  differential  equation  for  the 
matrix  T(t).  Since  the  elements  of  the  matrix  (dT(t)/dt)  are 
known.  the  elements  of  TCt)  are  integrated  numerically  from  the 
initial  values  specified  at  t=0. 

Siec  3:  With  the  elements  of  TCt]  tabulated  as  a  function  of 
time,  the  elements  of  U'(t]  are  calculated  as  a  function  of  time. 
The  matrix  U’(t)  =  T  U  T  (equation  85)  is  formed  for  a  given 
gamma- ray  frequency  . 

Step  i;  The  integral  of  equation  (89)  is  numerically  computed  as 

a  function  of  time,  the  magnitude  squared,  and  then  the  summation 

over  the  indices  of  the  excited  states  is  taken.  The  calculation 
®  2 

of  <  I  |c.Ct)|  >  as  a  function  of  time  is  taken  for 
k  =  3 

this  value  of  w  . 

Step  5:  This  summation  is  averaged  over  a  time  comparable  to 

several  natural  lifetimes. 

Step  S:  Steps  1  through  5  are  repeated  for  another  value  of 

until  the  spectrum  has  been  covered. 


A  program  which  performs  steps  1  through  3  has  been  written 


and  is  in  the  process  of  checkout.  Since  these  equations  are 
quite  complicated,  a  simple  case  was  chosen  which  could  also  be 
verified  without  resorting  to  numerical  analysis.  This  is  the 
static  case  where  M(t)  =  e^  or  ^  =  0°.  The  magnetization  lies 
in  the  xy  plane,  so  the  diagonal  elements  of  vanish. 

With  given,  and  the  splittings  ot  the  ground  and  excited 


levels  given  by 
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the  differential  equations  for  the  elements  of  T(t)  are 
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The  elements  of  the  T  matrix  in  this  special  case  are  found  to  be 
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Tgg  =  C-i/4j[y/3)  [sin  K^x  sin  3^x1  = 
where,  in  terms  of  the  parameters  of  the  computer  program. 
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All  other  T^  ■ 

=  0.  Tha 

t  T  is 

symmetric 

in  this 

case  is  due  to 

the  choice 

=  .M  =  M  . 

s 

This 

is  not  the 

usual  pr 

operty  of  T. 

The  matrices  TCt) 

and 

U’Ct)  have 

been 

calculated  and 

compared  with  the  functions  above.  There  is  very  good  agreement. 
Discrepancies  are  oscillatory  in  nature  and  co  not  seem  to  build 
up  in  time.  The  worst  case  deviations  of  U'Ct]  are  about  0.3%. 

Work  continues  on  this  approach  with  the  next  step 
being  the  checkout  of  each  computational  step. 
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SIGNIFICANCE 

To  dress  an  isomeric  state  requires  a  certain  arrangement  of 
nuclear  .levels  that  would  make  them  undetectable  to  conventional 
techniques  of  nuclear  spectroscopy.  Our  method  of  FMS  is  the 
only  means  found  to  date  that  can  be  used  to  search  for  this 
combination  among  the  29  best  candidates. 

In  order  to  advance  the  theory  of  FMS  during  the  current 
reporting  period  we  have  generated  a  computer  program  to  integ¬ 
rate  the  Landau-Li f shi tz  equation  for  the  foil  geometry.  Also  we 
are  in  the  process  of  completing  the  checkout  of  the  program 
described  in  section  5  which  theoretically  predicts  the  sideoand 
amplitudes  when  the  time  dependence  of  MCt),  given  by  the  solu¬ 
tion  of  the  Landau-Li f shi tz  equation,  is  known.  Since  this 
program  is  rather  general,  it  can  also  be  utilized  to  predict  the 
Mossbauer  absorption  spectrum  of  other  nuclear  species  embedded 
in  the  ferromagnetic  foil,  for  example  a  Sn  nucleus  embedded  in 
Fe.  Further  modifications  of  the  program  are  planned  to  extend 
its  capacity  so  that  more  general  motions  of  Mitj  may  be  incorpo¬ 
rated. 

The  successes  to  date  of  the  new  FMS  apparatus  for  nuclear 
spectroscopy  indicate  that  a  much  higher  resolution,  by  per.baps 
six  orders  of  magnitude,  can  be  achieved  through  a  reasonable 
upgrade  of  the  apparatus.  It  the  range  o!  tunability  does  extend 
to  the  ferromagnetic  spin  resonance  CFSRl  frequency,  then  it  will 
be  possible  to  construct  a  swept  frequency  device  capable  of 
continuously  tuning  over  a  range  of  10^^  linewidths,  an  enormous 
improvement  in  the  state-of-the-art  of  nuclear  spectroscopy. 
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Figure  1.  Tne  neart:  ol  tne  Frequency  Moduiatea  tipec  t  r  ome  c  a  r  irMo. 
IS  a  real  time  computer  intertace  with  a  linear  motor,  an  n i 
signal  generator,  ana  a  vray  detector.  Tne  v-rav  aoscroer, 
mounted  in  thin  cover  glass  slides  tor  rigidity,  is  suDjectec  to 
an  Hi  alternating  magnetic  tield.  'I'ne  computer  scans  tnrougn 
trequencies  ot  the  HI  tield  to  obtain  a  spectrum. 
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Figure  b.  Primed  coordinate  system 
detined  in  ttie  plane  ot  the  toil  sued 
chat  M  always  cooresponds  to  the  z  axis. 
The  vector  is  perpendicular  to  the 
plane  of  the  foil. 
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Figure  11.  The  Fourier  oecompos i c i on  ol 
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Figure  13A.  Coordinate  system 
tor  a  "rigid"  domain  wall 
motion  under  the  intluence  ot 
an  excernaiiy  applied  tield 


B  .  Adapted  from  Figure  2.4  ot 
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by  T .  H . 


O'Deii  . 


f-igure  13B.  Coordinate  system 
tor  domain  wail  motion  in 
which  the  magnetization  is 
tilted  slightly  out  ot  the  xz 
plane.  Analysis  snows  that 
the  moving  wall  sus.t  develop  a 
component  ot  M  which  lies 
along  the  direction  ot  mo¬ 
tion.  Adapted  trom  Figure  2.5 
ot  Fer  romaone  todvnami cs  by 
T.  H.  O’Dell^. 
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Figure  14.  The  magnetization  in  a  wail  region  o:  a  permaiioy 
tiim  0.1  urn  thick.  The  computation  ot  the  exact  behavior  ol 
domain  wails  in  conductors  is  based  on  the  minimization  ot  the 
interaction  energies,  including  the  magnetostatic  energy.  Ot 
note  is  the  tact  that  the  magnetization  always  lies  para.Liel  to 


the  surface  ot  the  conductor,  thus  minimizing  the  magnetostatic 


energy.  Adapted  from  A.  Aharoni 
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